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Abstract: We investigate Nj = 1 QCD in external magnetic fields on the lattice. The 
background field is introduced by means of the so-called Schrodinger functional. We adopt 
standard staggered fermions with constant bare mass am = 0.025 and magnetic helds with 
constant magnetic flux up to a^eH ~ 2.3562. We find that the the deconhnement and chiral 
symmetry restoration temperatures do not depend on the strength of the applied magnetic 
field. Our method allow us to easily study the effects of the external magnetic helds on the 
QCD thermodynamics. We determine the inhuences of applied magnetic helds to the free 
energy, pressure, and equation of state of strongly interacting matter. 
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1 Introduction 

Strong interactions are described by quantum chromodynamics (QCD), a local relativistic 
non-abelian quantum field theory which is not amenable to perturbation theory in the low- 
energy, large-distance regimes. However, many fundamental questions are linked to the 
large scale behavior of QCD. In particular, non-perturbative approaches to QCD can be 
used to account for the different phases of hadronic matter under extreme conditions. 
Recently, the study of the effects of strong magnetic fields on the QCD phase diagram 
has become a topic of increasing interest (for a recent review, see Refs. [1, 2]). In the non- 
perturbative regimes this problem can be efficiently approached by lattice QCD simulations 
with dynamical quarks. 

The study of lattice gauge theories with external background fields has been pioneered in 
Refs. [3, 4] for the U(l) Higgs model in an external electromagnetic field. In the continuum 
a background field can be introduced by writing: 

^ . ( 1 . 1 ) 

In the lattice approach one deals with link variables U^{x). Accordingly, on the lattice 
Eq. (1.1) becomes: 

U^{x) ^ U,{x)Ujr\x), ( 1 . 2 ) 

where U^^^{x) is the lattice version of the background field A^^{x). As a consequence the 
lattice action gets modified as: 

S[U] 5[[/] +6S[U,U^^^], (1.3) 

where 6 S[U, takes into account the influence of the external field. An alternative 
method, which is equivalent in the continuum limit, is based on the observation that an 
external background field can be introduced via an external current [5, 6]: 

’ ( 1 - 4 ) 
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so that the action gets modified as: 


S ^ S + Sb, 


(1.5) 


where: 

Sb = J dx J^^\x) A^{x) = - ^ J dx F^f{x)Fy^{x) . ( 1 . 6 ) 

The main disadvantage of these approaches resides on the lack of gauge invariance for non- 
abelian gauge theories. The issue of gauge invariance, however, does not pose if one is 
interested in QCD in external magnetic fields. In fact, let us consider the lattice partition 
function of QCD with / flavors of dynamical staggered quarks: 

z = j VU Yl[detM{U)]^ , (1.7) 

where Sg is the gauge field action and M is the fermion matrix for a staggered quark with 
bare mass amf. 

FIn,m,iU') — ^ ^ ^Ui/(^n)6jji^n+0 Uy(Tn'j6jyi,n—uj T OjTTl^ 5irYi,n t 

v=l 

V,{n) = (_i)-i+-+--i . (1.8) 


Since magnetic fields couple only to quarks, external magnetic fields can be introduced in 
the lattice action by replacing in the fermion mass matrix Eq. (1.8) the gauge field links 
according to Eq. (1.2), where the [/®^*( 3 :)’s are U(l) elements corresponding to the external 
magnetic fields with continuum gauge potential A^^^{x). For instance, if we consider con¬ 
stant magnetic fields directed along the X 3 direction, then the continuum gauge potential 
in the Landau gauge reads: 

Af^{x) = 4,2 xiH . (1.9) 

Therefore, we may write: 


t/“*(x) = C/|^*(x) = C/“*(x) = 1 , C/|^*(x) = cos{qfeHxi) -t- ism.{qfeHxi) , ( 1 - 10 ) 


where e is the (positive) elementary charge and ( 7 / is the quark charge (g^ = 2/3 , = 

— 1/3). Since the lattices have the topology of a torus, the magnetic field turns out to be 
quantized [4]: 


a^qf eH 


2tt 


^ext 5 


next integer 


( 1 . 11 ) 


where Lg is the lattice spatial size. Indeed, in the recent literature this approach has been 
adopted in extensive numerical simulations of QCD in external magnetic fields [7-23]. 

An alterative approach to put background fields on the lattice has been proposed since long 
time [24]. Indeed, that proposal allows to overcome the problem of gauge invariance by 
implementing background fields on the lattice by means of the manifestly gauge-invariant 
lattice Schrodinger functional. In this paper we present an exploratory study of lattice 
QCD immersed in a uniform external magnetic field. The background field is introduced 
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by using the Schrodinger functional. Moreover, for simplicity, we restrict ourself to one 
flavor of staggered dynamical quark. 

The plan of the paper is as follows. In Sect. 2, for completeness, we briefly discuss our 
method to introduce background fields on the lattice. In Sect. 3 we present the results of 
our numerical simulations for several local observables. We also address the problem of the 
possible dependence of the pseudoritical couplings on the magnetic field strengths. Sect. 4 is 
devoted to the discussion of the effects of magnetic fields on QCD thermodynamics. Finally, 
our conclusions are relegated in Sect. 5. 

2 Magnetic Fields within the Schrodinger Functional 

For reader’s convenience, let us briefly review background fields in lattice gauge theories 
within the Schrodinger functional. Firstly, we illustrate the method in pure gauge theories. 
In Ref. [24], to overcome the gauge invariance problem in presence of background fields, it 
was proposed that background fields on the lattice could be implemented by means of the 
gauge invariant lattice Schrodinger functional: 

= j ^ ( 2 . 1 ) 

where the functional integration is extended over links on a lattice with the hypertorus 
geometry and satisfying the constraints {xt = X 4 is the temporal coordinate) 

Uk{x)U=o = Ul^\x), (A: = 1,2,3). (2.2) 

One also imposes that links at the spatial boundaries are fixed according to Eq. (2.2). In 
fact, in the continuum this last condition amounts to the requirement that fluctuations over 
the background field vanish at infinity. This approach has been applied for both abelian 
and non-abelian gauge theories with different background fields [25-34]. 

The effects of dynamical fermions can be accounted for quite easily. Indeed, when includ¬ 
ing dynamical fermions, the lattice Schrodinger functional in presence of a static external 
background gauge field becomes: 

JUk (Lt ,x)=Uk (0,x)=i7“‘ (x) 

= [ VUe-^G det M , (2.3) 

JUk {Lt ,x)=Uk (0,x)=i7|"‘ (x) 

where Sp is the fermion action and M indicates the generic fermion matrix. Notice that 
the fermion fields are not constrained and the integration constraint is only relative to the 
gauge fields. This leads to the appearance of the gauge invariant fermion determinant after 
integration on the fermion fields. As usual we impose on fermion fields periodic boundary 
conditions in the spatial directions and anti-periodic boundary conditions in the temporal 
direction. In fact, Eq. (2.3) has been employed to study the dynamics of QCD with two 
degenerate staggered quarks [35, 36], as well as the quantum Hall effect in graphene [37]. 
In the case of QCD in constant magnetic fields the constraints in the lattice Schrodinger 
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functional need to be slightly modified to take into account that the magnetic held is 
coupled only to quarks. To this end, we impose that during the upgrade of the gauge links 
U^^{x) = I, while for the upgrade of the fermion helds U^^{x) = I X where: 

= 6k,2qfeHxi. (2.4) 

Since our Schrodinger functional Z[U^^^] is dehned on a lattice with periodic boundary 
conditions, usually we impose that: 

62{xi,X 2,X3,X4) = 92ixi + Ls,X2,X3,X4) . (2.5) 

As a consequence the magnetic held H turns out to be quantized: 

27r 

a^qf eH = — next (2.6) 

with next integer. However, it should be kept in mind that we are dealing with a periodic 
lattice with hxed boundary conditions, so that is is not strictly necessary to impose the 
quantization Eq. (2.6) and the "integer" next can be an arbitrary real number. 

3 Numerical Results 

We perform simulations of lattice QCD with one-havor of rooted staggered quark. Our 
numerical results were obtained by choosing as gauge action the Wilson action: 

Sg = ^ Re [Tr U^Ax)]) (3.1) 

where U^v{x) are the plaquettes in the (/r, n)-plane and (3 = \. Therefore we are led to 
consider the following lattice Schrodinger functional: 

Z[C/r*] = [ VU [detM(C/)]3 , (3.2) 

JUk{Lt,S)=Uk{0,x)=Ul^^(x) 

where the staggered fermion matrix is given by Eq. (1.8). To perform the functional integra¬ 
tion over the SU(3) links we have made use of the publicly available MILC code [38] which 
has been suitably modihed by us in order to introduce the boundary constraints Eq. (2.2). 
Ah simulations make use of the rational hybrid Monte Carlo (RHMC) algorithm. The func¬ 
tional integration is performed over the lattice links, but the links at the spatial boundaries 
are hxed according to Eq. (2.2). Accordingly, the links which are frozen are not evolved 
during the molecular dynamics trajectory and the corresponding conjugate momenta are 
set to zero. The length of each RHMC trajectory has been set to 1.0 in molecular dynamics 
time units. For each value of the gauge coupling f3 and the magnetic held eH we collected 
4000 - 5000 trajectories, and about 10000 trajectories around the critical coupling. To 
allow thermalization we typically discarded 1000 trajectories. The statistical errors were 
estimated by means of boostrap combined with binning. 

In the present exploratory study we consider lattices of size Lg = 24 and Lt = 4: and hxed 
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bare fermion mass mo = am = 0.025. At fixed Lt the temperature of the gauge system 
T = is changed by varying the coupling constant /3. 

Since the smallest quark electric charge is Igj = 1/3, from Eq. (2.6) we get: 

eH = ^ next (3.3) 

Different strengths of the external magnetic field are labelled by the parameter next ac¬ 
cording to Eq. (3.3). We performed simulations for next = 0,1,3, corresponding to field 
strength a?eH = 0, 0.7854, 2.3562 in lattice units, and assumed Qf = \ (up quark). Note 
that, the case of down quark Qf = —^ can be recovered with next = ~ ^ {a?eH = — 0.3927). 
In fact, to check the dependence of the free energy on the magnetic field we have also per¬ 
formed numerical simulations for next — ~ ^ ■ 

For the sake of completeness let us discuss, briefly, how the background magnetic 
field influences the dynamics of the gauge system. We said that to update the gauge 
system we used the rational hybrid Monte Carlo algorithm. As it is well known (see for 
instance Ref. [39]), to simulate the fermion determinant one introduces color-triplet scalar 
pseudofermion fields. The pseudofermion action depends on the inverse of the staggered 
fermion matrix Eq. (1.8). In the molecular dynamics one solves the equations of motion 
of the momenta conjugated to the gauge links. The derivative of the gauge momentum 
is called the force term, which is the formal derivative of the effective action with respect 
to the gauge potential. Thus, the force term consists of two contributions, namely the 
gauge force term and the fermion force term. Our boundary conditions correspond to set 
= j on the Xt = 0 hypersurface and at the spatial boundaries of the lattice in the 
gauge force term, while = I X in the fermion force term. To maintain the 

above constraints during the molecular dynamics, the momenta conjugated to the frozen 
gauge links are set to zero. 

3.1 Local Observables 

In this Section we are interested in the effects of the external magnetic field on several local 
observables. First, we consider the gauge action which, following the MILC convention, we 
define as: 

Gaction = ( X] Re Tr[/^j,(x)] \ . (3.4) 

In Fig. 1 we display the gauge action as a function on the gauge coupling (3 for three 
different values of the magnetic held. Since Gaction is a pure gauge quantity, it couples 
to the magnetic held only through quark loops. Therefore we expect that this quantity 
should manifest a very weak dependence on the magnetic held. Indeed, Fig. 1 shows that 
the effects of the magnetic held on the gauge action are at most of order 10“^ (see Fig. 5). 
Interestingly enough, we see that the gauge action increases as a function of eH in the 
strong coupling region whereas it decreases in the weak coupling region. In fact the three 
different curves displayed in Fig. 1 cross near the critical coupling /3c — 5.0. 

A more interesting quantity is the quark chiral condensate: 
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Figure 1. The gauge action Eq. (3.4) versus /3 for g/ = | and next =0,1,3. 




(3,5) 


which should display a pronunciate dependence on the magnetic field. In Fig. 2 we dis¬ 
play ('0 V’) versus the gauge coupling (5 for three different values of the magnetic field eH. 
The chiral condensate was computed by noise estimators with 4 random vectors. It is ev¬ 
ident that the chiral condensate increases as a function of eH for all temperatures. For 
comparison, in Fig. 2 we also display the real part of the Polyakov loop expectation value: 


(L) = 


1 


Lt 


3L3 


En Tr C/ 4 (x, xt) 

\ X xt=Q / 


(3.6) 


The Polyakov loop L, likewise the gauge action, is a pure gauge observable. Nevertheless, 
Fig. 2 shows that the Polyakov loop displays a sizable dependence on the magnetic field. In 
particular, we see that L increases with eH for all temperatures as for the chiral condensate. 
This behavior can be qualitatively understood if the quark free energy decreases with the 
applied magnetic field. In fact, later on we will show that the strongly interacting system 
behaves like a paramagnetic medium, i.e. positive magnetic susceptibility. 

Another interesting feature of Fig. 2 is the crossing of the chiral condensate and the Polyakov 
loop near the critical temperature. In fact. Fig. 2 seems to suggest that the pseudocritical 
gauge coupling /3c does not manifest a strong dependence on the magnetic field eH. 
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Figure 2. The chiral condensate Eq. (3.5) (left) and the real value of the Polyakov loop Eq. (3.6) 
(right) versus /3 for g/ = | and next =0,1,3. 


3.2 Pseudocritical couplings 


In this section we address the problem of the possible dependence of the pseudocritical 
coupling on the magnetic field. In general, the (pseudo)critical coupling is determined as 
the value for which some relevant susceptibilities exhibit a peak. In the present paper, to 
precisely localize the peak in the relevant susceptivity we parametrize the peak region with 
a Lorentzian function: 


m = 


ai 


(3.7) 


02 (/3 - /3c)^ + 1 

Our estimate of the critical coupling /3c is obtained by fitting the susceptivity to Eq. (3.7) in 
the peak region. We use the Polyakov loop susceptibility as well as the disconnected part of 
the chiral susceptibility to locate the transition temperature to the high temperature phase 
of QCD. 

First, we consider the disconnected chiral susceptibility: 


-i/; 


1 1 
Ll Lt 16 


(([Tr M-^f) 


[{Tr . 


(3.8) 


In Fig. 3 we show the disconnected chiral susceptibility as a function of the gauge coupling 
for three different values of eH. As usual the disconnected chiral susceptibility displays 
a sharp peak near the chiral critical coupling. Interestingly enough, the dependence of the 
chiral susceptibility on the magnetic field is almost relegated to the peak region. Moreover, 
the peak values increase with eH signaling that the chiral transition sharpens in presence 
of a non-zero magnetic held. Notwithstanding, we hnd that, within our statistical uncer¬ 
tainties, the chiral critical coupling does not depend on the magnetic held. 
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Figure 3. The disconnected chiral susceptibility Eq. (3.8) versus (3 for 9/ = | and next =0,1,3. 
The continuous lines are the results of the fits of the chiral susceptibilities to Eq. (3.7). The estimate 
critical couplings are reported in the legend. 

We have also considered the Polyakov loop susceptibility: 

Xi = Ll {{L^) - (Lf) . (3.9) 

Results for the Polyakov loop susceptibility are shown in Fig. 4. In this case we see that 
the dependence of the Polyakov loop susceptibility on the magnetic field is less pronounced 
with respect to the disconnected chiral susceptibility. Moreover, Fig. 4 shows that the peak 
of XL decreases with increasing eH. This means that the applied magnetic field tends to 
smooth out the deconfinement transition. However, even in this case the deconfinement 
critical couplings does not display any appreciable dependence on eH. Moreover, we find 
that the chiral and deconfinemnet critical couplings agree for any values of the magnetic 
field strengths considered in this work. Finally, as further check, we have considered the 
variation of the gauge action: 

action = Gaction{LH ^ 0) — Gaction{LH = 0) , (3.10) 

which is known to display a peak in the critical region. Of course, by using AGacUon we 
may estimate the pseudocritical couplings only for non-zero magnetic field strengths. 

In Fig. 5 we show AGacUon versus /3 for the three different values of the magnetic field 
employed in the present work. As discussed in Sec. 3.1 the effects of the magnetic field on 
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Figure 4. The Polyakov loop susceptibility Eq. (3.9) versus (3 for Q/ = | and next = 0,1,3. The 
continuous lines are the results of the fits of the Polyakov loop susceptibilities to Eq. (3.7). The 
estimate critical couplings are reported in the legend. 


Table 1. Summary of the values of the critical couplings /3c estimated with different operators for 
the magnetic field strenghts considered in this work. 


Operator 

a?eH 

/3c 

^disc 

0 

5.047(5) 

0.7854 

5.044(6) 


2.3562 

5.049(4) 


-0.3927 

5.046(2) 

XL 

0 

5.060(5) 


0.7854 

5.048(10) 


2.3562 

5.041(11) 


-0.3927 

5.051(5) 

^Gaction 

0 

- 


0.7854 

5.078(6) 


2.3562 

5.079(5) 


-0.3927 

5.070(16) 


the gauge action are tiny. Moreover, the non monotonic dependence of the gauge action 
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Figure 5. AGacUon Eq. (3.9) versus /3 for g/ = | next =1,3, and g/ = -1 next = 1- The 
continuous lines are the results of the fits of AGacUon to Eq. (3.7). The estimate critical couplings 
are reported in the legend. 


on eH is clearly displayed in Fig. 5. In any case, we see that AGacUon does display a 
well developed peak in the critical region. We find that the peaks in AGacUon ai'e located 
at a systematically slightly larger values of the gauge coupling with respect to the chiral 
and Polyakov loop susceptibilities. In a hnite volume this is, of course, not unexpected. 
Indeed, we recall that we are using lattices with hxed boundary conditions and, as previous 
studies showed, the gauge action turns out to be more susceptible to hnite volume effects. 
Nonetheless, what it is relevant is that the pseudocritical couplings do not depend on the 
magnetic held eH. 

For reader convenience, in Table 1 we summarize our estimates of the critical couplings /3c 
as function of the magnetic held eH. From Table 1 we may safely conclude that the critical 
temperature does not depend on the external magnetic held. 


4 Thermodynamics in external Magnetic Fields 

The partition function Z\U^^] allows us to dehne observables that can be used to establish 
the equation of state of the theory. Such observables play an important role in describing 
the thermodynamic properties of the system. 
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The free energy density is related to the logarithm of the partition function as: 


f{T,H) = - I logZpn = log^[C/r] • (4.1) 

The pressure is given by the derivative of TlogZ\U^^\ with respect to the volume. As¬ 
suming that we have a large, homogeneous system, differentiation with respect to V is 
equivalent to dividing by the volume. Therefore in the thermodynamic limit the pressure 
can be written as minus the free energy density: 

p{T,H) = - f{T,H) . (4.2) 


Using the well-known relation between the trace anomaly (also called interaction measure) 
and the derivative of the pressure: 



= e{T,H) - 3p{T,H) = ^ 

p{T,H) 

2^4 ’ 

(4.3) 

one can easily calculate the energy density: 




= I{T,H) + 3p{T,H) , 


(4.4) 

the entropy density: 

siT,H)+ p{T,H) 
s[T,H) = ^ 


(4.5) 

and the speed of sound: 

2 _ dp 
^ de ■ 


(4.6) 


As usual, we need to renormalize the free energy density by subtracting the divergent 
zero-point energy. To do this it is enough to subtract the zero temperature contribution. 
Thus we dehne: 


fr{T,H) = f{T,H) - f{0,H) , Pr{T,H) = -fr{T,H). (4.7) 

The zero temperature contributions are conventionally obtained by performing simulations 
on lattices with Lf = Lg. Moreover, since we are interested in the thermal magnetic 
properties of our system, we will focus on: 

Afr{T,H) = fr{T,H) - fr{T,H = 0) , Apr{T,H) = -Afr{T,H). (4.8) 

In a Monte Carlo simulation, one cannot compute the partition function directly. The most 
frequently used method in practice is the integral method, in which a derivative of the free 
energy with respect to some parameter serves as observable, which then gets integrated 
again to yield the free energy density. Since we are doing simulations at hxed Lt , it is 
convenient to take derivatives with respect to the bare gauge coupling /3. The expectation 
values of the derivatives with respect to /3 of our partition function correspond to the average 
Wilson action. Thus, we have: 

7 ^ -^ - j^{{{Sw)h,t - {Sw)h,o) - {{Sw)h=o,t - {Sw)h=o,o)} ■ (4.9) 
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Figure 6. Eq. (4.9), versus T/Tc for g/ = | next =1,3, and qf = -1 next = 1- 


In Fig. 6 we report our results for the /3-derivative of Apr{T, H) (normalized to T^) versus 
the ratio T/Tc for three different values of the magnetic field. We recall that the temperature 
corresponding to a given value of the gauge coupling is given by the relation T = 

For the dependence of the lattice spacing on the gauge coupling we used the two-loop 
/3-function. Accordingly, we have: 


«(/5) Aqcd = /qcd(/3) , 


(4.10) 


where /qcd(/ 3) is the asymptotic scaling function of QCD with one dynamical fermion 
Nf = l-. 


/qcd(^) = ( 


-bi/{2bl) 


exp 




11 - iNf 

, bo = 3 / 


12bo J ’ (47r)2 ’ ^ (47r)^ 

For definitiveness, the critical temperatures have been obtained by using the pseudocritical 
gauge coupling estimated by means of the chiral susceptibilities. Obviously, a direct deter¬ 
mination of the physical scale should be preferable. However, it is known that the lattice 
violations to the asymptotic scaling law Eq. (4.10) are within a few percent. So that, the 
adopted approximation is adequate to the purpose of the present exploratory study. 

From the derivative of the pressure, after numerical integration, we may easily obtain 
Apr{T, H) = —Afr{T,H). In Fig. 7 we display the normalized pressure Apr(T, H) ver¬ 
sus the temperature for three different values of the magnetic field. Fig. 7 shows that the 


102 - f W 

3 ^ . (4.11) 


12 - 














Figure 7. versus T/Tc for g/ = § Uext =1,3, and qf = —^ Uext = 1- 


magnetic contributions to the renormalized pressure is clearly different from zero even for 
T < Tc, and it seems to vanish rapidly for low temperatures. This behavior can be naturally 
accounted for within the Hadron Resonance Gas model (see, for instance, Ref. [40, 41]). 
On the other hand, for temperatures above the critical temperature the pressure increases 
rapidly in qualitative agreements with perturbative calculations in the high-temperature 
regime [42]. 

As concern the energy density, using Eq. (4.3) we may write: 

I{T,H) ^ Aer{T,H) - 3Apr{T,H) ^ ^ _d_ 

T4 dT d(5 


Apr{T,H) 

2^4 


(4.12) 


dT 



Rm , 


(4.13) 


where Rp{(d) is the lattice /Jgc'D-function. According to our approximation, we have: 


Rm 


/qcd(/3) 

^/qcd(/3) 


(4.14) 


Using Eqs. (4.12), (4.13), and (4.14) we determined the so-called interaction measure dis¬ 
played in Fig. 8 for three different values of the magnetic field. 
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Figure 8. The interaction measure , Eq. (4.12), versus T/Tc for qf = ^ next =1,3, and 

Qf — 3 ^ext — 1 ■ 


After that, the magnetic contributions to the renormalized energy density can be 
straightforwardly obtained as: 


Asr {T,H) AI{T,H) 




J14 


+ 3 


Apr{T,H) 

rp4 


(4.15) 


In Fig. 9 we show the renormalized energy density. Even for the magnetic contribution to 
the energy density we find two different regimes for T < (conhned phase) and T > 
(deconfined phase). Having determined the magnetic contribution to the pressure and 
energy density, in principle one can construct the equation of state and obtain the entropy 
density and the speed of sound by means of Eqs. (4.5) and (4.6). For the purposes of 
the present paper we do not discuss any further this matter. We, merely, observe that 
in the deconfined phase T > Tc the magnetic contribute to the energy density increases 
slower with respect to the pressure by increasing the temperature. This behavior leads to 
a stiffening of the equation of state. 

Let us, finally, address the problem of the magnetic susceptibility. As is well known, for 
small magnetic field strengths we may write for the free energy density (see, for instance. 
Ref. [43]): 

A fr{T, H) = - ^ Xmag{T) , (4.16) 
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Figure 9. , Eq. (4.15), versus T/T^ for qf = ^ Uext =1,3, and qf = — ^ next = 1- 


where Xmag is the magnetic susceptibility. Therefore, to determine the magnetic suscep¬ 
tibility we need to check if A fr{T, H) = —Apr{T,H) scales with at least for small 
enough magnetic field strengths. To this end, in Fig 10 we display for different 

values of the magnetic field strength. In fact, we see that the free energy density seems 
to scale with for {qf = next = 1)) and {qf = —next = 1)) and for temperatures 
not too far from the critical temperature. We note, however, that in the high-temperature 
region physical observables are more affected by finite volume and cutoff effects. On the 
other hand, we see clearly that for the strongest magnetic field used in this paper {qf = |, 
n-ext = 3) the scaling with is badly violated. It is useful to give the corresponding 
values of the magnetic field in physical units. To this end, we use the known flavor depen¬ 
dence of the QCD critical temperature reported in Ref. [44] to infer that for Nf = 1 the 
critical temperature is Tc ~ 190 MeV. This corresponds to a lattice spacing a ~ 0.26fm. 
So that, in the critical region, we estimate y/\eH\ ~ 0.48, 0.67, l.lTGeV corresponding to 
next = — 0.5 , 1, 3 respectively. Thus, we see that for magnetic field strengths not exceeding 
1.0 GeR the free energy density seems to display an approximate scaling with within 
our statistical uncertainties. This is in qualitative agreement with the results in Ref. [7] 
where the strongest magnetic field used was y|e77| ~ 0.85GeV. In any case, for physical 
applications, we recall that the magnetic fields relevant for heavy-ion collision experiments 
are of order ^J\eH\ ~ 0.1 GeV. 
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Figure 10. versus T/T^ for g/ = -^ next = 1, and g/ = f next =1,3. 

To determine the magnetic susceptibility we are lead to consider the free energy density for 
\/\eH\ < 1.0 GeV where we can safely apply Eq. (4.16). Indeed, in Fig. 11 we report our de¬ 
termination of the magnetic susceptibility as function of the temperature. For comparison, 
we also display our determination of the magnetic susceptibility for ^/\eH\ ~ 1.17 GeV. 
From Fig. 11 we see that the magnetic susceptibility is positive in the whole temperature 
range explored in the present study. Moreover, the magnetic susceptibility increases mono- 
tonically with the temperature. Therefore, the strongly interacting medium behaves as a 
paramagnetic substance both below and above the critical temperature Tc- It is remark¬ 
able that our results for the magnetic susceptibility is in fair qualitative and quantitative 
agreement with Ref. [19] where it has been considered Nf = 2 + 1 QGD with physical 
quark masses, discretized on a lattice by stout improved staggered fermions and a tree level 
improved Symanzik pure gauge action. 

5 Conclusions and Discussion 

In conclusion, let us summarize briefly the main results of the present paper. We investi¬ 
gated QGD with one flavor of staggered quark in an external magnetic field on the lattice. 
The external magnetic field has been introduced by means of the so-called Schrodinger 
functional. We have investigated the magnetic properties of one-flavour quarks and gluons 
in thermal equilibrium for magnetic field strengths up to ■\J\eH\ < 1.17 GeV. In particu- 
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Figure 11. The magnetic susceptibility as a function of the temperature for different 


magnetic field strengths. 


lar, we focused on the effects of the magnetic field on several local observables and found 
results in qualitative agreement with recent results in the literature obtained with a dif¬ 
ferent method, as described in Sect. 1, to implement external magnetic fields in QCD on 
the lattice. We have, also, looked for the effects of the magnetic field on the critical tem¬ 
perature. Surprisingly, we found that the critical temperature does not change even for 
the strongest magnetic field used in the present work. This is in striking contrast with the 
results in the literature. However, since we used one flavor rooted staggered quark which is 
known to be strongly affected by taste symmetry violation effects, one could suspect that 
our results on the critical temperature is merely due to lattice artifacts. Indeed, presently 
we are simulating the same physical system by adopting highly improved staggered quarks 
(HISQ) where the taste symmetry violations are dramatically reduced. Nevertheless, our 
preliminary simulations adopting HISQ quarks to do not yet display a clear dependence of 
the pseudocritical temperature on the background magnetic field. In any event, we plan to 
report progress on this subject in a future paper. 

We evaluated the magnetic contributions to the pressure, energy density, and free energy. 
Our results are in qualitative agreement with previous investigations. In particular, we con¬ 
firm that the free energy density scales with for small enough magnetic field strengths. 
Moreover, we determined the magnetic susceptibility and found that the strongly inter¬ 
acting medium behaves like a paramagnetic substance both below and above the critical 
temperature in agreement with previous results in the literature. 
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